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Abstract 



Explicit expressions for optical tomograms of the photon-added coherent states, even/odd photon- 
added coherent states and photon-added thermal states are given in terms of Hermite polynomials. 
Suggestions for experimental homodync detection of the considered photon states are presented. 
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Different kinds of nonclassical photon states were studied during last two decades (see, e.g. reviews 
[1]). Among the nonclassical states there are so called photon-added and photon subtracted states [2j. 
These states were considered in [3], [3], [5]. In [5] it was shown that in photon subtracted states the 
y—i mean number of photons may be either less or greater than of initial state depending on the photon 
statistics, while in photon-added states the meal number of photons is always greater then that of initial 
state. Experimentally these states were studied in [7J, [8], [9], [10], [11] where the commutation relation 
of bosonic annigilation and creation operators was considered. The experiments were done by using the 
homodyne detection procedure of photon states applied in [12] , [13] to reconstruct the Wigner function 
of photon states. On the other hand the optical tomograms of the quantum state can be considered as 
primary object which contain complete information of the quantum states and can serve as alternative of 
density matrix |14] , [15] . In view of this the homodyne detection of photon states used in experiments 
with photon-added states is interesting per se. The homodyne detection provides as experimental output 
the optical tomogram which is fair probability distribution of homodyne quadrature X depending on extra 
parameter 9. This angle parameter is so called local oscillator phase which is controlled parameter varying 
in domain > 9 > 2tt. Note, that this probability satisfies the symmetry relation w(X, 9+ir) = w(—X, 6), 
and thus in the range > 9 > ir it completely describes the quantum state. 

In [7] the authors used conditional stimulated parametric down-conversion in a nonlinear optical 
crystal for the production of the photon-added states. Here one high-energy pump photon can annihilate 
into two photons that obey the global energy and momentum conservation laws and are normally emitted 
into symmetrically oriented directions (signal and idler modes). To generate photon-added state one has 
to inject a seed (initial) state into a signal mode of the parametric amplifier. The conditional preparation 
of the target photon added state will be take place every time that a single photon is detected in the 
correlated idler mode. 

The primary light source for the experiment is a mode-locked laser whose pulses are frequency-doubled 
to become the pump for degenerate parametric down-conversion in a type-I beta-barium borate crystal. 

The Wigner functions of single photon added coherent states were reconstructed for experimentally 
measured data for |a| = 0.03, 0.9, 2.6. 
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In [S] the authors have probed quantum commutation rules by addition and subtraction of single 
photons to/from a light field. They used the same procedure to generate photon-added, photon sub- 
tracted, photon-added-then subtracted, and photon-subtracted-then-added states, which are different 
from each other and from the original seed thermal state. They experimentally measured the quadrature 
distribution function of photon-added and photon-subtracted states for initial thermal state with mean 
number of photons n = 0.57. Note, that this measured distribution function coincides with optical to- 
mogram because the density matrix of this state exhibits no phase dependence. In development of this 
work in [10] the authors measured tomograms of superpositions of photon-added-then-subtracted and 
photon-subtracted-then-added states , depending on phase between two these states. 

In [9] the parametric down-conversion in a nonlinear crystal also used to generate photon-added 
coherent states, and non-Gaussianity of these states were investigated. Experimental Wigner functions 
were reconstructed for single photon coherent states with \a\ = 0, 0.48, 1.02. 

In [11] Experimental preparation of single photon added thermal states and reconstruction of the 
Glauber-Sudarshan quasiprobability (P-function) for these states have been done. 

It has been shown |16j that the higher order photon-added coherent states \a,m) corresponding to 
m > 1, also can be realized in the parametric down-conversion scheme. 

Generation of the photon-added coherent states also is possible in the cavity-atom interaction [2]. 

More over, the interaction in optomechanical systems (micro-resonator interacting with laser light) 
by proper detuning can be used to generate photon-added coherent states in the optomechanical domain 

m- 

A common feature of the mentioned experimental techniques is that they are all used interacting 
bipartite systems. During evolution, the two subsystems are entangled. This enables to make suitable 
conditional measurements on one of the subsystems so that the other subsystem is prepared in a photon- 
added state. 

In report [E] a comparative study of two processes, namely, parametric down-conversion and atom- 
cavity interaction, that can generate photon-added states is presented. By expressing the time-evolved 
states in a suitable non-orthogonal basis, it is established that the later method generates ideal photon- 
added states, a feature that is not present in the atom-cavity scheme. Further, the parametric down- 
conversion itself is shown to be capable of generating ideal m-photon-added state, without requiring 
higher order processes. 

The aim of this work is to study optical tomograms of different kinds of photon-added nonclassical 
states. We consider photon-added coherent states and thermal states. Even/odd coherent states [19] 
with added photons will be also studied in this paper. The obtained below optical tomograms could be 
used to compare the experimental results of the homodyne detection of the photon states, i.e. the output 
experimental optical tomogram (e.g. in [7J, [H], [H]) [ID]) with the tomograms calculated in this our paper. 
It is worth to note that some recent studies of photon-added states are presented in [18J. 

We also want to point out that the consideration of homodyne detection of one-mode photon-added 
states in the tomographic probability representation and the analysis of the experimental precision of 
the obtained optical tomograms [7j, [8], [9], [TO], [H] has the aspects which have to be discussed in the 
connection with other experiments where different photon states were studied, see, e.g. |20j . 

In all the published results not all the tests related to properties of the optical tomograms associated 
with checking the uncertainty relations and studying the dependence of optical tomogram on local oscil- 
lator phase have been given in detail and in explicit form. Another goal of our work is to suggest to pay 
more attention to the variety of the possible tests which permit to do better estimation of accuracy of 
the measuring such primary object as optical tomogram identified with quantum state. 

The paper is organized as follows In Sec. 2 short review of optical tomographic representation of 
the photon states is given. In Sec. 3 we review a model of parametric oscillator and photon-added states 
appeared in this model are considered. The resume of results is presented in conclusion Sec. 4 
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2 Optical tomographic representation 



In this section we review the tomographic probability representation of quantum states |15j . Given a 
density operator p (one mode photon state). The symplectic tomogram M(X, p,, v) which is normalized 
probability distribution of homodyne quadrature X depending on two real parameters p and u is defind 



M(X,p,u) = Tx{p5(X - pq-vp)} = {X,p,u\p\X,p,u), 



(1) 



where q and p are operators of the first and second photon quadratures, and \X, p,v) is an eigenvector 
of the hermitian operator pq — vp for the eigenvalue X. For arbitrary state with wave function ^>(y) the 
tomogram is expressed in terms of fractional Fourier transform of the wave function 



M(X,p,u) = (X,ti,v\y){*\X,fi,v) 



2ttW 



ip 2 iXy 

— y 



(2) 



In experiments with homodyne detecting photon states the optical tomogram is measured. This 
tomogram reads 

w(X,9) = Tr{p5(X - qcos9 - psin9)} = (X,6\p\X,0), (3) 



where \X, 9) is an eigenvector of the hermitian operator qcos9 — psin9 for the eigenvalue X. 
The tomograms 

@ and([3]) are related 



w{X,0) = M(X, cos 6, sin ( 



M(X,p,v) 



X 



w 



\Jp 2 + v 2 



, tan 



p 



The density operator can be reconstructed as [21] 



— / M{X,p,v) j( x -w-P») dX dp dv. 
2vr J 



(4) 



(5) 



(6) 



In view of §5§ the density operators can be determined in terms of the optical tomogram w(X, 9). The 
invertible relations of the measurable tomograms and density operator provide possibility to use the 
optical tomogram as primary object containing complete information of quantum states. In fact, the 
quantum evolution equation and energy level equation were written for optical tomogram in explicit 
form in |22j . The quadrature statistics can be obtained from optical tomogram 



(X n }(9) = J X n w(X,9) dX. 



Also the photon number statistics can be obtained from optical tomogram. In fact, 



(q n ) = J X n w(X,9 = 0) dX, 



(p n ) = J X n w{X,9 = vr/2) dX. 

Also mean photon number reads 

{n) = \j X 2 [w(X, 9 = 0)+ w(X, 9 = vr/2)] dX - l -. 



(7) 



(8) 



(9) 



(10) 
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The experimental tomograms must satisfy uncertainty relations. These relations in the tomographic 
form read 

J X 2 w(X, 6 = 0) dX - (^j X w(X, 9 = 0) dX^j 

J X 2 w(X, 9 = tt/2) dX - (J X w(X, 6 = n/2) dx\ 



111) 



For example for coherent state the product in this relation equals 1/4. 

We point out that the optical tomography can provide possibility of direct measuring the purity of 
photon quantum state [23] , 



3 Optical tomography of photon-added states 

3.1 Photon-added coherent states 

Let us find tomograms of studied in [3j time dependent photon- added coherent states \a,m,t) of one 
mode parametric oscillator with the Hamiltonian 

h = ^+ n 2 (*)f , m = i. 

The state \a, m, t) is defined as follows: 

\a,m,t}= U(t)\a,m,0) = U(t)\a,m) = (m!L m (-|a| 2 ))~ 1/2 U(t)a +m \a), (12) 

where L m (z) = Lm \z) is the Laguerre polynomial [2H[25], \a) is the initial coherent state, U(t) is the 
unitary evolution operator 

U(t) U + (t) = 1, 17(0) = i. (13) 
As shown in [3] the expression for the state \a,m,t) in the coordinate representation has the form 

.* \ m/2 / 



(q\a,m,t) = (m\L m (-\a\ )) 



-1/2 / e 



a (q\a,t), 



2eJ y\e\ V 2e 

where H m (z) is the Hermite polynomial [241 125j . (q\a,t) is the time-dependent coherent state 

( q \ a , t) = exp + - S£ - H!) , 

which was considered, e.g. in [26J, c-number function e{t) satisfies the equation 

s(t) + n\t) e(t) = 0, 
with the initial conditions e(0) = 1, e(0) = i, which means that the Wronskian is 

EE* — E*E = —2%. 



(14) 



(15) 



(16) 



(17) 



The state \a, m, t) is pure, thus we can find tomogram of it from the formular Q. After some calculations 
we get 

-1 f / ^- , \ „ . i /o ^ 2 



M am (X,/j,,v,t) 



(m\L m (-\a\ 2 )Y 

^7T2 m \EU + 



exp 



\a\' 
~2 



X 1 



Xe + iV2av fe* \ ( |e| 2 (//e + ue) 
\e\(n£ + vi) V 2e a J\ e 2 (fie* + ve*) 

^ 2 



1/2' 



+ 



V2aX a 2 e* 



+ 



iva 



2\/j,e + ve\ 2 fie + ue 2e e(fie + ve 



(18) 



4 



The substitutions fi = cos 9 and v = sin 9 to (18) gives us according to Q the optical tomogram 
w am (X,9,t). In the case of stationary Hamiltonian (fl(t) = 1, e = e lt ) we have 



w am (X,9,t,Q = 1) 



(mLM-M 2 ))- 1 



vr2 r 



Hm ( X 



a 

7=2 



-i(t+6) 



exp |-X 2 - \a\ 2 + 2v / 2XRe(ae 



-i(t+6i)> 



Re(a 2 e- 2j ^)} 



(19) 




3.2 Photon-added even/odd coherent states 

Photon-added even/odd coherent states \a±,m) were introduced in [4J as follows 

a jm \a±) . . 

({a±\a m ai m \a±)) 1 -< z 

where \a±) is an even/odd coherent states |19j . |27j and m is an integer. The normalized states \a±) can 
be written as 

e M 2 /2 

\ a +) = —7= . | , 9 (l«) + I -a)), 
za/cosIi |a| z 

e M 2 /2 

|a_) = — (|a) — \ —a)). 

2ysinh|a| 2 

Photon-added even/odd coherent states \a±,m) in terms of photon-added coherent states \a,m) are 
represented as follows [3J: 

I \ ( eH 2 Lm (-|q| 2 ) V /2 

|a± ' m)= 2(eH 2 L m (-|a| 2 )±e-H 2 L m (|a| 2 )) (|a,m> ± | - a,m» (21) 
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This formula in view of ([I]) and (12) enables us to get the result 

M <x±m (X,[i,i/,t) ■ 



eH 2 L m (-|a| 



2(eH 2 L m (-|a|2)± e -H 2 L m (|a| 2 )) 
x {M am (X,ix,i/,t) + M- am (X,fj,,u,t) 
± 2Re [{X, fj,, u\a, m, t){— a, m, t\X, /J,, u)]} , 



(22) 



and with respect to (|4j) we have w a±m (X,6,t) = M a±m (X, fi = cos9,v = sin6>,t), and in the case of 
stationary Hamiltonian with f2(t) = 1 we have w a±m (X, 8, t) = w a±m (X, t + 6). 




a) 




b) 



Figure 2: 

a = 0.1 (a] 



Optical tomograms w a+ (X, 9) of even photon-added coherent states for Q{t) = 1, m = 1, 
, a = 1 (b). 



3.3 Photon-added thermal states 

A thermal state is the most classical state of light, formed by a statistical mixture of coherent states. 
The density operator of thermal state has the form px = Z~ 1 e~ H l T , Z = Tve~ H / T . The associated 

optical tomogram reads w^(X,9) = — = e~ x / 2<J , where a 2 = ^coth^. The density matrix of the 
photon-added thermal states in representation of the Fock states \n,t) is given by 

(k,t\p T rn\n,t) = (k,t\^ m p T a m \n,t) = 5 k J -> — e -(«" m )/ T , (k, n>m). (23) 

m! [n — my. 

The symplectic tomogram of this state can be written as follows 

oo 

M Tm {X,n,v,t) = ^J(n + m,t\p T m\n + m,t)\(X, n,v\n + m,t)\ 2 . (24) 

n=0 
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a) b) 

Figure 3: Optical tomograms w a _(X, 9) of odd photon-added coherent states for O(t) = 1, m 
a = 0.1 (a), a = 1 (b). 



Noting that \(X, /j,, u\n, t)\ = M a= o }n (X, [i, u, t), with the help of (18) we can find the result 

/i e -l/T\m+l 

l i_e ) „-X 2 /\n£+ve\ 2 



M Tm (X,^,u,t) = — — — — ■ e 

y/TT ml 2 m |/x£ + vc 



~ e -n/T 
X > — , 

^ n! 2 n 



n+m 



Xe ( \e\ 2 (fie + z^e) x ' 



n=0 

For the stationary Hamiltonian the optical tomogram can be read 



e\(fj,£ + ve) \£ 2 (fj>£* + ue 



(25) 



= ^m! 2™ 6 2^ ^ • W 



n=0 



Thus the optical tomogram of the photon-added thermal state for the stationary Hamiltonian is not 
depend on time and on the parameter 0. This formula with the help of the known summing expression 
for Hermite polynomials [21] 

n _ „2wi/2 PYn ( 2xy V -(x 2 + y 2 ) V 2 \ _ ^ H n {x)H n (y) 
U V ) exp I I - ^ , r? , 

v ' 7 n=0 

is transformed to the differential expression for WT m (X,6) = WTm(X) 
It gives us for m = 1 and m = 2 the following relations: 

WTm=1 {x) ~ ^/T^T- w Wwf + r^f f (28) 
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U>Tm=2(X) 



(1 - e- 1 ^)- 
[ 4X 4 



-X 2 tanh(l/2T) 



-2/T 



+ 



4X 2 (2e" 1 / T - 1) 



+ 



2e" 2 / T + 1 



|(l + e- 1 /T)4 ( 1 + e -l/T)2( 1 _ e -2/T) (!_ e -2/T)2 

In Figures [T] - [4] we illustrated tomograms of photon-added states for some examples. 



(29) 
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a) 



b) 



Figure 4: Optical tomograms WT m {X, 6) of thermal states for T = 1, m = 1 (a), m = 2 (b). 



4 Conclusion 

To resume we point out the main results of our paper. We calculated the optical tomograms of 
photon-added states. The initial states to which the photon are added are chosen as coherent state, 
thermal state, even and odd coherent state. We considered these states since they have been realized 
experimentally in [7j,[8], |10j . We pointed out that the optical tomograms in probability representation 
of quantum mechanics can be considered as primary objects containing complete information on quantum 
states (see, e. g. [15J). 

In view of this the reconstructing the Wigner function or Husimi function do not add extra information 
on the state property. Thus the problem is to measure the optical tomogram with highest possible 
accuracy. The precision of the experiments for measuring optical tomogram can be checked using such 
criteria as quadrature uncertainty relation with purity (temperature) dependent bound [28J. This criteria 
was suggested in |23| to control precision of the homodyne experiments. In the mentioned experimental 
works [7] , [8] , [9] , [IU] , [H] , including the homodyne two- mode studies [20] the checking correspondence 
to the above criteria has not been done. It is especially important, taking into account that the local 
oscillator phase contribution can influence the final results if the states under the study differ from the 
states (thermal or photon-added thermal) with optical tomograms not depending on the local oscillator 
phase. 

It is of special interest to test the experimental optical tomogram in the sense that it must satisfy 
the discussed in introduction symmetry condition with respect to shift of local oscillator phase by n. 
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In view of this we suggest to measure accurately the optical tomograms not only in domain of local 
oscillator phase > > tt but to investigate accurately also the domain it > 8 > 2tt to check the 
symmetry properties of the tomogram. Using the theoretically calculated optical homodyne probabilities 
as a standard for comparison with the experimental results, one can estimate measure of necessity of the 
precision improvement of the technique in future experiments. 
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